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ABSTRACT 

We investigate the color-magnitude relation for globular clusters (GCs) — the so-called "blue tilt" 
— detected in the ACS Fornax Cluster Survey and using the combined sample of GCs from the ACS 
Fornax and Virgo Cluster Surveys. We find a tilt of = d{g — z)/dz — — 0.0257 ± 0.0050 for the 
full GC sample of the Fornax Cluster Survey («5800 GCs). This is slightly shallower than the value 
7^ = -0.0459 ± 0.0048 found for the Virgo Cluster Survey GC sample (willOO GCs). The slope for 
the merged Fornax and Virgo datasets (wlGOOO GCs) is 7z = —0.0293 ± 0.0085, corresponding to a 
mass-metallicity relation of Z cx jV^o.43±o.i2^ g^^^ ihdX the blue tilt sets in at masses in excess of 
M 2 X lO^M©. The tilt is stronger for GCs belonging to high-mass galaxies (M* > 5 x lO^^M©) than 
for those in low-mass galaxies (M, < 5 x IO^^Mq). It is also more pronounced for GCs with smaller 
galactocentric distances. Our findings suggest a range of mass-metallicity relations Zqq oc A/q^^" '' 
which vary as a function of host galaxy mass/luminosity, a scaling similar to that observed for dwarf 
spheroidal galaxies. We compare our observations to a recent model of star cluster self-enrichment 
with generally favorable results. We suggest that, within the context of this model, the proto-cluster 
clouds out of which the GCs formed may have had density profiles slightly steeper than isothermal 
and/or star formation efficiencies somewhat below 0.3. We caution, however, that the significantly 
different appearance of the CMDs defined by the GC systems associated with galaxies of similar mass 
and morphological type pose a challenge to any single mechanism, including self-enrichment, that 
seeks to explain generically the observed GC color-magnitude relations. We therefore suggest that 
the detailed (and stochastic) merger/accretion histories of individual galaxies have likely played a 
non-negligible role determining the distribution of GCs in the color-magnitude diagrams of individual 
GC systems. 

Subject headings: galaxies: clusters: individual: Fornax - galaxies: clusters: individual: Virgo - 
galaxies: fundamental parameters - globular clusters: general 
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1. INTRODUCTION 

In recent years, deep space-based imaging, in com- 
bination with high-resolution spectroscopy of individ- 
ual stars, has revealed the existence of multiple stel- 
lar populations in several globular clusters (GC_s] be- 
longing to t he Mi lky Way (e.g., Lee e t i 
Bedin et al' (120041) . P iotto et al. ^2m$_A 
Villanova et al. ()2007D . Milone et al. (|200"8[ )). This 
finding has both undermined the long-held view that 
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star clusters are simple systems consisting of single- 
age, single-metallicity stellar populations, and renewed 
interest i n GC self-enrichment scenarios (e.g., Fran k 
& Gisler (fl976l ). Smith Plel), Gnedin et al. 
Parmentier & Gil more ( ^001), Dopita & S mith 0986), 
Morgan fc La ke |l989} ). Thoul et al (I2002D . Par- 
mentier (|2004D Recchi & Danzi ger (I2005D . Caloi & 
D'Antona (|2007t ). Bailin & Harris (|2009r )). 

As a complement to such studies, observations of (un- 
resolved) GCs in galaxies beyond the Local Group of- 
fer the advantage of much larger and more homogeneous 
databases of GC colors and magnitudes with which to 
investigate overall trends between integrated GC prop- 
erties and the properties of the GC systems as a whole. 
Observations from HST, using both WFPC2 and ACS , 
have been pi votal in this area (e.g.. Cote et al. (2004), 
Harris eioLdaOOeD & ([2009ai) . Peng et al. (2006), Jordan 
et al. ()2007a| )). One of the more surprising results to have 
emerged from these studies has been the discovery of a re- 
lationship between the magnitudes and colors of individ- 
ual GCs, in the sense that GCs associated with the blue 
subpopulation becom e prog ressively redder at high lumi- 
nosities (Harris et al. ()2006l) , S trade r et al. (l2006f ). Mieske 
et al. ()2006[ ). Spitler ei al. ()2006[ ). Humph rey (|2009D . 
Cockcroft et al. (|2009| ). Harris et al. (|2009b[ ). Forbes et 
aL (HpiS))- 

This relation, which is now commonly refered to as 
the "blue tilt", has also been detected in ground-based 
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imaging of extr agalac tic GC sys tems (Forte et al. ()2007l ). 
Wehner et al. (|2008l ). Harris (|2009D\ and possibly for 



M31 as well (Fan, Ma & Zhou ([20091)). If interpreted 
as a mass-metallicity relation, the tilt corresponds to a 
scaling with mass of Z oc ^Vf*''^""-^. This is compara- 
ble to the relation Z oc defined by the faintest 
dwarf sp heroid al galaxies of the M ilky Way (e.g. Simon 
& Geha (|2007t) , Kirby ei a/. ()2008D ) which have luminous 
masses comparable to GCs but gravitating masses 1-2 or - 
der of magnitude lar ger. Recently, Mieske et al. (|2006l) 
and Blakeslee et al. ( 2010[ ) have shown that a unimodal 
GC metallicity distribution with an underlying mass- 
metallicity trend may create a color magnitude rela- 
tion with a distinct blue color peak — a consequence 
of the non-linear co lor-me tallicity relation for GCs (see 
also Richtler et al. (|2005[ ). Yoon et al. (2006), Peng et 
al. (120061) '). 

The most homogeneous and extensive sample of extra- 
galactic GCs currently available is that from t he AC S 
Virgo Cluster Survey (ACS VCS, Cote et al. (jlOOl)). 
From imaging of 100 Virgo cluster early-type galaxies, 
more than 1 0000 GC candidates were i dentified (e.g., 
Peng etoL^M), Jordan et al. (|2009( )). In Mieske 
et al. (|2006l ). ACSVCS XIV, we examined the color- 
magnitude relations of these GCs. A highly signifi- 
cant correlation, 7^ = d(g-z)/dz = -0.0374 ± 0.004 
{Z oc M° ''*), between color and magnitude was found 
for the subpopulation of blue GCs in the co-added sam- 
ples of the three brightest Virgo cluster galaxies (M49, 
M87 and M60). In general, the trend was found to be 
more pronounced for GCs belonging to more massive 
host galaxies, but there were also clear galaxy-to-galaxy 
differences; significant blue tilts were found for the GC 
systems of M87 an d M60 but not for th at of M49 (see 
also Strader et al. ()2006D ). Peng et al. ( |2009f ) have re- 
cently confirmed the presence of a blue tilt in M87 using 
much deeper (50 orbit) ACS imaging in the F606 {V) 
and F814 (/) filters, in contradiction to a claim based 
on these same data th at the blue tilt is an observational 
artifact (Waters et al. (j2009| )). 

In this paper, we analyze the color-magnitude relation 
of GCs in the ACS F ornax Cluster Survey (ACS FCS, 
Jordan et al. (|2007a[ )). This survey, which imaged 43 
early-type galaxies belonging to the Fornax cluster, is 
nearly identical in design to the ACS VCS, although it 
targeted galaxies in a new and different environment. We 
also analyze the combined GC samples from the Fornax 
and Virgo surveys, to improve the overall statistics and 
highlight possible cluster-to-cluster differences. In what 
follows, we focus on four main topics: (1) the variation in 
the tilt slope as a function of environment; (2) the mass 
scale where the tilt first appears; (3) the constraints im- 
posed on GC self-enrichment models posed by our find- 
ings; and (4) the amount of tilt that may arise from GC 
dynamical evolution. 

2. SELECTION OF GLOBULAR CLUSTERS FROM THE 
ACS FORNAX CLUSTER SURVEY 

The ACS FCS sample consists of 43 early-type (E, SO, 
dE, dE,N, dSO) members of the Fornax cluster. Each 
galaxy was imaged in the F475W and F850LP filters 
(~ ff475 and zsso) for a total of 760 and 1220 sec, respec- 
tively. This filter combination gives roughly a factor- 
of-two improvement in wavelength baseline and metal- 



licity sensitivity compared to the "canonical" {V — I) 
color index (Cote et al. 2004). The identification of 
bona-fide GCs from these images is performed in the 
size-magnitude plane as described in Peng et al. (|2006[ ) 
and Jordan et al. (2009). This selection procedure — 
which is possible because the half-lig ht rad ii of GCs are 
marginally resolved (Jordan et al. ()2005[ )) at t he di s- 
tance of Fornax (d = 20.0 Mpc; Blakeslee et al. ()2009[ )) 
— greatly reduces contamination from both foreground 
stars and background galaxies. 

The ACS images have been reduced using a ded i- 
cated p ipelin e that i s described in Jordan et al. ( ()2004a| ). 
(|2004bD and (|2007a[ ). In brief, the reductions consist of 
image combination, galaxy modeling, model subtraction, 
rejection of obvious background galaxies, and the mea- 
surement of magnitudes and sizes for candidate GCs u s- 
ing the program KINGPHOT (Jordan et al. (|2005D ). 
The result is a catalog containing integrated g and z 
magnitudes, (g — z) colors and half-light radii, r^, for 
all candidate GCs. Magnitudes and colors are corrected 
for foreground extinction using the reddening maps of 
Schlegel et al. (199^ and ex tinctio n ratios for a G2 star 
as specified in Sirianni et al (|2005[) . 

To estimate the contamination by background galax- 
ies, an identical reduction procedure was applied to 
17 blank, high-latitude fields observed by ACS in the 
F475W and F850LP filters (see Peng et al. (|2006[ ) and 
Jordan et al. ((2009i ) for details). The relative distri- 
butions of contaminants and GCs in the size-magnitude 
plane then make it possible to assign a GC probabil- 
ity, "Pgc, to every source detected in each field. In our 
analysis, all sources with Vgc > 0.5 are considered to 
be GCs. The residual contamination from background 
galaxies that are assigne d "Pgc > 0.5 is negligible for this 
study (see Mieske et al. ||20M ). ACSVCS XIV in the fol- 
lowing) . 

3. RESULTS 

Fig. [T] shows color-magnitude diagrams (CMDs) for 
GCs in the Fornax cluster in vs. {g — z): for NGC 
1399 alone, the entire FCS sample, and the FCS sample 
separated by host galaxy magnitude at Mb = —19.75 
mag. This magnitude corresponds to a stellar mass of 
Af, « 5 X 10^°Mq (Peng et al. 2008). Although some- 
what arbitrary, this choice divides the GC samples into 
two nearly equal halves, while at the same time corre- 
sponding roughly to several noteworthy mass scales in 
the physics of galaxy formation: i.e., the approximate 
transition between "young and old" galaxies, at ~ 
3 X IO^^Mq (e.g., Kauffmann et al. 2003, 2004), and 
the division between galaxies showing central luminosity 
deficits and excesses in their surface brightness profiles, 
at - (5-10)x101"Mq (Ferrarese et al. 2006ab; Cote 
et al. 2006, 2007). Throughout this paper, we will re- 
fer to the galaxy sample with Mb < —19.75 mag as 
the high-mass sample, and that with AIb > —19.75 as 
low-mass sample. 

We use the SBF dis tance moduU from Mei et al. (j2007D 
and Blakeslee et al. (|2009f ). which have an internal pre- 
cision of « 0.08 mag, to assign absolute magnitudes to 
the associated GCs. Fig. [2] shows the CMDs of the full 
VCS sample, the combined FCS and VCS sample, and 
the combined FCS and VCS sample after separation into 
high-mass and low-mass galaxies. 
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3.1. KMM Fitting Procedure 

Anal ogous to ACSVCS XIV, we apply KMM (Ashman 
et al. (|1994D ) in heteroscedastic mode to the CMDs of 
these sub-samples, subdivided into 25 luminosity bins 
containing the same number of GCs. By adopting a 
fixed number of luminosity bins, we are able to get a 
straightforward comparison of the color scatter between 
different subsamples. The only exception is the CMD 
of NGC 1399 in Fig. [J for which we adopt 10 luminos- 
ity bins due to the smaller number of GCs. For each 
bin, KMM is used to determine the best-fit blue and red 
peak magnitudes of a bimodal Gaussian color distribu- 
tion. In the heteroscedastic mode, we allow for Gaussians 
of different widths to be fitted to the blue and red GC 
subpopulations. 

To obtain an estimate of the uncertainty involved in 
the choice of the initial peak position guesses, we perform 
the KMM fitting in two iterations. We first run KMM using 
the mean peak colors of the full samples (i.e., FCS, VCS, 
or FCS and VCS combined) as initial guesses. From 
the KMM results for a given subsample, we then calcu- 
late the peak color {g — z)o averaged over all luminosity 
bins, and calculate the dispersion, do, of the peak col- 
ors around their mean. Typically, ctq was found to be 
in the range 0.05 to 0.1 mag. For the second iteration, 
we then run KMM twice, with two different initial guesses: 
(g — z)o ± (Tq. For the final measurement we adopt the 
mean of the KMM results obtained from the two initial 
guesses. In Figs. [Hand [2] these fitted positions of the blue 
and red peaks are plotted over the respective CMDs. The 
mean is indicated by a filled circle, while the results from 
the two initial guesses are shown as asterisks. In most 
cases, the KMM results were found to be identical for the 
two initial guesses. We also note that for all the sam- 
ples investigated in this paper, bimodal Gaussian color 
distributions are strongly prefered over unimodal Gaus- 
sian distributions (confidence levels >99%). Note that in 
ACSVCS XIV, it was shown that the KMM results are re- 
liable for Mz-{g — z), while for Mg-{g — z) they are biased 
toward progressively bluer colors at bright luminosities 
(due to the depopulation of the red peak at high g-band 
luminosities). We therefore do not consider Mg-{g — z) 
CMDs. 

In these plots, we show as vertical dashed lines the 
mean color expected for the blue and red peak based on 
the relation between host galaxy magnitude and GC peak 
colors found by bi nning independent color distributions 
(Peng et al. (|2006D ). The mean colors of the blue peaks 
agree very well with the KMM fit results in this paper. For 
the full VCS, and VCS plus FCS, samples, the KMM fits 
and median colors give a small redward offset of < 0.05 
ma g for t he red peak compared to the results of Peng et 
al. ( 20061 ). due to the slightly asymm etric s hape of this 
peak (see Figs. 2 and 5 of Peng et al. ()2006D ). 

To measure the amplitude of the tilt, we perform least- 
square fits of the KMM data points for the blue and red 
peaks as a function of magnitude, excluding GC candi- 
dates brighter than Mz — —12.5 mag and magnitude 
bins fainter than Mz = —8.1 mag. The resulting lin- 
ear relations are indicated in the CMDs of Figs. [1] and 
[21 The errors on the fit are calculated by resampling 
the points using as the dispersion the observed scatter 
around the fitted relation. The measured slopes and er- 



rors ^z = d{g ~ z)/dz are shown in Table [T] We also in- 
dicate an estimate of the corresponding mass-metallicity 
scaling Z cx M". This is obtained from converting the 
(g — z) datapoints to [Fe/H] using a biquadratic transfor- 
mation based on Peng et al. (2006) and fitting a linear 
relation to the data points in magnitude- [Fe/H] space. 

The slope of this transformation is '^^gl^j — 5 dex/mag 

in the colour range of the blue peak {{g — z) ~ 0.95 mag) 

and — 2 dex/mag in the colour range of the red 

peak {{g — z) ~ 1.35 mag). We assume a constant mass- 
to-light ratio over the luminosity range investigated. 

As a check on the KMM fits, we calculate the slope using 
the median colors blueward and redward of the limiting 
color between the blue and red peak. In Figs. [T] and [2l 
the median values are indicated as open triangles, and 
the fit to them as a dotted line. Note that we define the 
limiting color as the point where the contribution from 
the blue and red peak is equal; this is generally different 
from the minimum of the color distribution. The limiting 
color is first determined for each luminosity bin based on 
the results from KMM. We then fit a linear relation to the 
limiting color as a function of magnitude and use this 
limiting color (the long dashed lines in Figs. [T] and [H to 
separate the blue and red peaks. We report the slopes of 
these median estimates in Table [H These generally agree 
quite well with the KMM results. Only the median in the 
most luminous bin is somewhat bluer than the KMM re- 
sult. This is because the KMM fit in the brightest bin typ- 
ically shows a redward "jump" of the color distribution 
with respect to fainter magnitudes (see Sect. [?75| due to 
the near merging of the blue and red peak, an effect that 
is not as apparent using the median estimate due to the 
smoothly changing colour limit adopted for the median 
calculation. The blue tilt measured from the median 
is hence, on average, 8-10% weaker than that measured 
from KMM. For the CMDs of NGC 1399 and the Fornax 
high-mass sample — which exhibit the most dramatic 
tilts — we also indicate the "minimum" color-magnitude 
relations obtained from determining the median colors 
blueward and redward of a magnitude-independent lim- 
iting color. 

3.2. Measured Slopes for Various Subsamples 

Table [T] shows that the slope obtained for NGC 1399 
is quite dramatic, with jz = -0.0878 ± 0.0250. This 
trend is significant at roughly the 3.5tT level. In mass- 
metallicity space, this corresponds to a scaling relation 
of Z oc MO-82±o.23 rpj^g morphology of the CMD of 
NGC 1399 shows that this strong slope originates from 
a near-merging of the blue and red pea ks at the high- 
est luminosities (see also Dirsch et al. (|2OO30 : Bassino 
et al. (j20061) ). Nevertheless, we emphasize that KMM still 
gives a confidence level of 98 — 99% for prefering a bi- 
modal Gaussian over a unimodal Gaussian for the bright- 
est bins. 

When considering the full sample of GCs from the 
FCS, we find a shallower but still highly significant tilt 
of jz = -0.0257 ± 0.0050. This is weaker than the tih 
jz = -0.0459 ± 0.0048 found for the fuU VCS data set 
(see Fig. [2] and Table [TJ. For the subsample of GCs 
belonging to the high-mass FCS galaxies, we find a 
slope of "fz = -0.0697 ± 0.0116. GCs belonging to the 
low-mass FCS galaxies show a rather weaker slope of 



4 



Mieske et al. 



7^ ^ -0.0186 ± 0.0066. This result for the FCS con- 
firms the finding from ACSVCS XIV that the slope be- 
comes weaker for GCs belonging to less luminous galax- 
ies. This also holds when considering the alternative 
fit results based on median colors, and even when us- 
ing the "minimum" slope derived in the previous subsec- 
tion. To reinforce this point, the bottom left panel of 
Fig. [2] shows the combined FCS and VCS sample, sub- 
divided into high-mass and low-mass galaxies. The tilt 
for the high-mass sample is 7^ = —0.0566 ± 0.0076, 
almost three times steeper than that for the low-mass 
sample, 7^ = —0.0210 ± 0.0034. The two slopes are in- 
consistent with each other at the 4.3cr level, and imply 
a range of mass metallicity scalings Z oc M^-^~^-'^ as a 
function of galaxy luminosity/mass. For the full sample 
of GCs from the two surveys (Fig. [2]), we find a slope of 
7^ = —0.0293 ±0.0085. This slope, which is, as expected, 
intermediate to those found for the FCS and VCS sam- 
ples, corresponds to Z oc jVf 0-43±o.i2^ 

The bottom right panel in Fig.[2]shows two GC subsets 
in the high-mass galaxies from the combined FCS and 
VCS samples, separated at a projected galactocentric dis- 
tance d = 5.5 kpc (corresponding to w 0.25-li?e for these 
galaxies; Ferrarese et al. 2006a, 2010). The slope found 
for the inner sample is 7^ = —0.0750 ± 0.0120, signifi- 
cantly higher than the value = — 0.0428±0.0071 found 
for the outer sample. The tilt difference is even more pro- 
nounced when using the alternative median fitting. This 
confirms the findings from ACSVCS XIV that the tilt is 
most significant for GCs lying at small projected galac- 
tocentric distances. 

3.2.1. M^P Revisited: Blue Tilt Differences Between 
High- Luminosity Galaxies 

Both in ACSVCS XIV and in Strader et al. (|2006l ) it 
was noted that no blue tilt was detected for the brightest 
Virgo cluster galaxy M49. To put this intriguing result 
into perspective, we compare in Fig. [3] the CMDs for 
NGC 1399 (Fornax), M49 (Virgo), and M87 (Virgo). We 
choose a subdivision into 15 bins for the KMM fitting 
to aU three galaxies. Both NGC 1399 and M87 show 
significant blue t ilts (A CSVCS XIV, Strader et al. (2006), 
and Peng et al. ()2009[ )). while M49 apparently shows no 
obvious color-magnitude relation. Comparing the overall 
morphology of the CMDs, it is clear that in the case of 
M49, the lack of bright, intermediate-color GCs is at least 
partly responsible for the lack of a tilt: i.e., the blue peak 
in M49 is narrow and there is a clear gap between the 
blue and red peaks over the entire magnitude range. 

This difference relative to M87 and NGC 1399 is 
difficult to reconcile within a scenario where a GC 
mass-metallicity relation is created entirely by a single 
mechanism operating within individual GCs (e.g., self- 
enrichment; see below). In all likelihood, the detailed 
accretion and merger histories of individual galaxies will 
also play a role in establishing the distribution of GCs 
within the color-magnitude plane. 

3.3. A Luminosity /Mass Threshold for the Blue Tilt 

It has been suggested that the blue tilt is driven mainly 
by GCs that are brighter than t he GCLF turnover mag- 
nit ude, (e .g. Harris et al. (j2006t ). ACSVCS XIV, Peng et 
al. ()2009t )). To investigate the luminosity /mass thresh- 
old of the blue tilt, we show in the upper left panel of 



Fig. |4] the dependence of slope 7z on the lower magni- 
tude limit used in the fitting procedure. In doing so, 
we have fitted a linear relation to the KMM peak results 
for the combined Fornax and Virgo samples, subdivided 
into high-mass and low-mass host galaxies (bottom left 
panel of Fig. [2]). The faint magnitude cutoff has been 
varied from Mz.cutjamt = -7.5 mag to -9.5 mag. 

For the high-mass sample, there is a notable trend for 
the derived slope to increase steadily from 7^ ~ —0.045 
to —0.075 as the faint cutoff magnitude is varied from 
cut, faint = —7.5 mag to —9.5 mag. This may be a 
signature of a non-linear color-magnitude relation at the 
upper end of the GC luminosity function. The open cir- 
cles in this plot show the change in slope when the bright- 
est magnitude bin, at = —11.2 mag, is excluded. In 
this case, the slope is roughly constant at 7^ ~ —0.045, a 
value of that is largely independent of the faint cutoff for 
Mz^cut, faint ^ —8.0. Hence, the steepening slope seems 
to be due to the increasing weight of the very brightest 
GCs, with < -11 mag. 

We conclude from this exercise that the blue tilt is 
firmly in place above M^^cut ~ —8 mag (1.5 — 2 x lO^M©) 
and further steepens at < —11 mag (> 2 — 3 x 
IO^Mq). Nevertheless, the upper left panel of Fig. H] 
shows clearly that the strong blue tilt for GCs in the 
high-mass sample is not driven by the brightest GCs 
alone since the open circles are consistently above the 
crosses. 

This plot also demonstrates once again that the blue 
tilt is more pronounced in the Virgo data. In this con- 
text it is crucial to note that the GCs in the VCS sample 
originate, on average, from more luminous galaxies than 
in the FCS sample, see the upper right panel of Fig. SI 
there is a difference of 0.5-1.0 mag in host galaxy magni- 
tude. The GC number ratios between the high-mass 
and low-mass samples is 43:57 for Fornax, but 60:40 
for Virgo. Given the overall stronger tilt found for the 
high-mass samples, this difference in host galaxy distri- 
bution appears to be the root cause of the stronger tilt 
found using the Virgo sample. 

In the bottom panels of Fig. |4] we investigate in de- 
tail how "fz changes as a function of the bright cutoff 
magnitude. We consider the high-mass and low-mass 
subsamples in the combined FCS and VCS sample and 
keep Mz = —8.1 mag fixed as the faint cutoff. In or- 
der to obtain a finer grid of peak positions at bright 
luminosities we: (1) apply a varying bright cutoff to 
the individual GC data points, instead of a cut to the 
KMM peak positions as above for the faint magnitude cut- 
off^^; and (2) adopt a slightly smaller bin size of 200 GCs. 
We then vary the upper cutoff in steps of 0.25 mag be- 
tween Mz = —12.5 mag an Mz = —8.5 mag. For each 
restricted sample, we fit red and blue KMM peaks. In 
the bottom left panel, we overplot the KMM peak results 
for all adopted bright GC magnitude limits. For the 
high-mass sample, the peak positions of the blue GC 
population "bend" slightly at higher GC masses. In the 
bottom right panel of Fig. IH this behavior is quantified 
by plotting 7^ as a function of the upper cutoff magni- 
tude. For GCs belonging in the high-mass sample, the 
slope is to within the errors independent of the upper 

The latter is adequate due to the much denser spacing of 
magnitude bins at fainter magnitudes 
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cutoff for —11 < Afz^cut, bright < —9 mag, but increases 
somewhat from 72 — —0.040 to — —0.052 when in- 
cluding also the brightest GCs with -12.5 < < -11 
mag. For GCs in the low-mass sample, no steepening 
for bright GCs is apparent. 

In summary, the tilt becomes stronger for the most 
massive GCs (> 2 - 3 x IO^Mq) belonging to the 
high-mass sample, but not for the GCs in the low-mass 
sample. 

3.4. Evidence for a Red Tilt? 

We now examine the evidence for a "red tilt" in the 
data: i.e., a color- magnitude relation for the subpopula- 
tion of red GCs. Table[T]shows that there is indeed a 2-3ct 
effect that is detectable in the high-mass samples, indi- 
cating a gentle mass metallicity relation of Z oc M"-^""-^. 
There is, however, no significant red tilt detected in the 
full sample of all GCs. Given this and the fact that the 
colour distribution in the red peak i s gene rally asymmet- 
ric (see Figs. 2 and 5 of Peng et al. (2006)), the evidence 
for a red tilt should still be considered preliminary. We 
note that a weak tilt is broadly consistent with the GC 
self-enrichment scenario, where the pre-enrichment level 
of red (metal-rich) GCs is significantly higher than that 
of the blue (metal-poor) GCs, thereby making the detec- 
tion of any self-enrichm ent m ore difficult (e.g. ACSVC S 
XIV, Strader & Smith (l2008h . Bailin & Harris ^OM)). 

We conclude this section with a note of caution regarding 
the detection of tilts, red or blue: the colors of individ- 
ual GCs typically scatter significantly about the individ- 
ual peak positions by 0.1-0.2 mag. This scatter - likely 
caused by stochastic variation in the GC pre-enrichment 
level - will allow the detection of any mass-metallicity 
trend only for systems of at least a few hundred GCs. 
The significance of the blue tilt detection in the GC sy s- 
tem of the disk galaxy NGC 5170 (Forbes et al. (|2010l) ) 
- which has about 600 GCs - is ^3a, showing that it 
would be difficult to detect subtle tilts in poorer GC 
system s like that of the Milky Way (see also Strader & 
Smith (|2008( ) and Baihn & Harris (2009.) ). 

4. DISCUSSION 

In ACSVCS XIV, it was shown that neither stochastic 
fluctuations within the stellar CMDs, contamination of 
the GC samples by tidally-stripped nuclei, nor the ac- 
cretion of GCs like those found in present-day galaxies 
could fully explain the observed color-magnit ude r ela- 
tions. Subsequently, in Mieske & Baumgardt ()2007D . it 
was shown through N-body simulations that field star 
capture is much too inefficient to account for any tilt. In 
this section, we examine the s elf-en richment scenario put 
forward by Bailin & Harris ()2009l ). including new con- 
straints imposed on this scenario by the observed trends 
with environment (§4.1). We then quantify the effects of 
GC dynamical evolution on the distribution of clusters 
within the color-magnitude plane (§4.2). 

4.1. Constraints on Self- Enrichment Scenarios 

Models for the self-enrichment of star cluster s hav e 
a long histo ry (see, e.g., Frank fc Gisler (|1976l ). 
Smith (|1996l) Gn edin et al. (|2002D Par mentier 
& Gilmore (1200 iD. Dopita & Smith (I1986D . Mor- 
gan fc Lake (|1989l ). Thoul et al. (|2002D . Recchi 



& Danziger (|2005D . Parmentier ()2004D . Baumgardt et 
al. ()2008l )) and they appear to offer a promising explana- 
tion for the blue tilt. The mass-metallicity relations as 
reported in this paper range from Z cx M°-^ to oc M°-^. 
This overlaps with the Z oc jVf0.6-o.7 relations fo und fo r 
dwarf sphero idal g alaxies (e.g., Simon fc Geha (|200l . 
Kirby et al. (|2008| )) which have luminosities, metallici- 
ties, and central velocity dispersions that are similar to 
those of GCs (although the latter appear to be baryon, 
rather than dark matt er, do minated). 

In Strader fc Smit h (|2008D . Bailin fc Harris (|2009t ) and 
Harris et al. ()2009b[ ) the self-enrichment scenario of star 
clusters is discussed with a focus on the blue tilt. The 
underlying assumption is that self-enrichment by SNII 
ejecta starts to become efficient when the energy output 
by SNII is comparable to the binding energy of the clus- 
ter's primordial gas cloud. This implicitly assumes that 
100% of the SN ejecta energy is converted to k inetic en- 
ergy of the cloud gas (see Tenorio-Tag le et a l. ()2007D for 
an alternative view). Bailin & Harris ()2009() parameter- 
ize the ability for self-enrichment of a proto-cluster cloud 
based on the metal-retention fraction, fz (equation 7 of 
their paper), 

logZc/Zrs, = log(10^'"-' + (1) 

where /, denotes the star formation efficiency. The par- 
ticular value of 0.38 used in the exponent of the sec- 
ond term is determined by the SN yield. For the metal- 
retention fraction, Bailin & Harris derive (equation 28 of 
their paper): 



fz = cxp(- 



EsNf*n 

102 MqGMc 



(2) 



This expression holds for the assumption of a 3D mass 
density profile of the protocluster cloud that falls off as 
p oc r~2 (i.e., an isothermal sphere), with rt being the 
truncation radius. In this equation, Esn is the typical 
energy released by a single SNII (= 10^^ erg) and Mc is 
the mass of the proto-cluster cloud. A Salpeter-like stel- 
lar mass function at the high-mass end is assumed. [Note 
that the present-day cluster mass is linked to the proto- 
cluster cloud mass via the relation Mc.today — Mcf*] 
For a more general density profile p on r^^ with /3 ^ 2, 
it holds for fz (equation 36 of their paper): 



fz = [1 - 



EsNf*rt{2 



(3) 



102 MqGMc 

Equations [5] and [3] show that the primordial star for- 
mation efficiency, /*, the proto-cluster mass-radius re- 
lation, and density profile, p, will determine the self- 
enrichment capability of a given cluster. Our aim is 
to constrain these parameters based on the observed 
GC mass-metallicity relation. To this end, we plot 
in the top panels of Fig. [5] GC metallicity vs. abso- 
lute magnitude and mass, separated into high-mass and 
low-mass samples. The [Z/H] values are obtained di- 
rectly from converting the (17 — z) datapoints to [Fe/H] 
using the b iquadratic transformation based on Peng 
et al. ()2006[ ). We also show separately the data for the 
three high-luminosity galaxies NGC 1399, M49 and M87 
(see §3.2.1). 

Fig. [5] shows the expected behaviour of log Zc/Zq as 
a function of mass for various choices of model param- 
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eters. The short dashed curve corresponds to th e ref- 
erence model adop ted by Baihn & Harris (|2009( ) (see 
also Harris et al. ()2009b[ )). For this model, the den- 
sity profile is taken to be an isothermal sphere with a 
mass-independent truncation radius of = 1 pc, and 
a star formation efficiency of f» = 0.3 (e.g. Boily & 
Kroupa (|2003D . Lada & Lada (|2003[ ). Baumgardt et al. 
2008). 

This particular set of model parameters does not repro- 
duce the data for either the high-mass nor the low-mass 
sample. Most notably, the observations point to an on- 
set of the blue tilt at w 2 x 10^ A/©, while the refer- 
ence model predicts self-enrichment to become impor- 
tant only above « 2 x IO^AIq. Moreover, the relation 
between mass and metallicity predicted by the model 
in the self-enrichment mass regime is Z oc 
This is significantly steeper than the observed relation 
of Z oc Af 0-3-o.7_ 

Thus, the data require a lower mass for the onset of 
self-enrichment than in the reference model, but at the 
same time, a weaker scaling between metallicity and mass 
in the self-enrichment regime. One way of achieving this 
is to assume a radial density profile that is steeper than 
an isothermal sphere: Fig. [S] shows that the particular 
choice oi (3 — 2.35 matches very well the observed mass- 
metallicity scaling for the high-mass sample. For such 
a profile, self-enrichment is able to start at lower cloud 
masses due to the deeper potential well in the cloud cen- 
ter, but its efficiency increases more gradually with in- 
creasing cluster mass (given that only a small fraction 
of gas is located close to the edge of the proto-cluster 
cloud). For the low-mass sample, the data require a vari- 
ation of P with cloud mass in the sense that lower-mass 
clouds have a steeper profile than higher-mass clusters; 
the shallower tilt for this sample can otherwise not be 
reproduced. The solid curve in Fig. [5] corresponds to the 
particular choice oi (3 = 2.3 -I- 0.23[6.2 — log(mass)]. 

An alternative way to match the observations for the 
high -mass sample is to adopt a lower star formation ef- 
ficiency of ~ 0.15. This implies that the progenitor 
cloud for a GC of a given mass was twice as massive as 
in the reference model, yielding a correspondingly larger 
metal retention fraction, and hence, a lower present-day 
onset mass for self-enrichment. At the same time, the 
choice of a lower reduces the overall efficiency of self 
enrichment, that is, the slope between [Z/H] and mass 
(Equ. [I}. For the low-mass sample, however, the curves 
with lower star formation efficiency do not reproduce the 
observed scaling of metallicity with mass. 

It is important to bear in mind that the assumption 
of a constant cloud radius is not necessary to reproduce 
the observations. Fig. [S] shows that clouds with con- 
stant density (i.e., r oc AI^^^) can also have levels of 
self-enrichment that are in agreement with the observa- 
tions. 

Finally, we note that radiative feedback from mas- 
sive stars also contributes to regulating star formation, 
with the total radiative energy input from OB stars 
into the ISM being more than an order of magnitude 
low er tha n the energy injected by SNII (Baumgardt et 
al. ()200 8f)'). Despite the lower amount of energy involved, 
radiative feedback may be an important factor for halt- 
ing star formation in massive clusters (e.g. Krumholz 



& Matzner ()2009[ ). Fall et al. ()2010[ )) in the few million 
years before the first SNII explode. This effect is not 
considered in the self-enrichment scenario discussed in 
the present paper, since we focus on the reprocessing of 
metals ejected by SNII, which dominate the budget of 
metal-enriched material in a star cluster. 

4.1.1. Environmental Dependence 

We have shown that the mass-metallicity scaling rela- 
tion, and the difference between GCs in high-mass and 
low-mass samples, may be explained if the density pro- 
files of the proto-cluster clouds are generally steeper than 
isothermal spheres, and if the density profiles for the 
proto-cluster clouds in the low-mass sample depended 
weakly on mass. Another effect that may contribute to 
the observed trends is that, for larger cluster masses, an 
increasing fraction of proto-cluster clouds with low star 
format ion ef ficiency survive their formation (Baumgardt 
et al. (|2008[ )'). This is possible since, for higher-mass 
clusters, gas expulsion will tend to happen more adia- 
batically and the tidal radius is larger, so that clusters 
will have more room for expansion folUowing gas expul- 
sion. This possibility is also supported by the finding of 
Baumgardt et al. ([2008) that, in order to turn a power- 
law cluster mass function into the co mmonly observered 
log- normal one (e.g., Jordan et al. (|2007bf ) and refer- 
ences therein) , a larger fraction of high mass clusters has 
to survive gas expulsion. Environmental variations in 
the importance of these effects are certainly conceivable: 
e.g., in the sense that metal retainment is favoured in the 
central parts of high-mass galaxies with external pres- 
sure confinement, different from regions in lower mass 
gala xies w ith a weaker external pressure (e.g., Parmen- 
tier (|2004[ ). Recchi & Danziger ([2005)). These variations 
may well infiuence the resulting mass-metallicity relation 
of GCs. 

However, such smooth trends can hardly explain the 
gross differences between the GC CMD morphology for 
galaxies of similar morphological type and luminosity, 
as is the case for M49, M87 and NGC 1399 (Fig. [3] 
and Fig. [5]). Instead, these differences may indicate 
that individual merger or accretion eve nts ha ve altered 
their G C pop ulations (e.g.. Cote et al. (|1998D or Hilker 
et al. (|1999( )). and hence contributed to modulating 
the blue tilt. This modulation might take the form of 
intermediate-age GCs formed in gas-rich mergers, the 
dissipationless accretion of GCs belonging to ga s-poor 
galaxies, or both (see the re views of West et al. (l200i ) 
and Brodie & Strader ()2006[ )). 

An alternative explanation for the variations in the 
color-magnitude relations among different subsamples is 
considered by Blakeslee et al. (|2010f l. These authors show 
that the non-linear relationship between GC color and 
metallicity may lead to different tilts in color space be- 
tween samples of different mean metallicity, even though 
the underlying GC mass-metallicity relations may be 
quite similar. In the present paper we adopt a color- 
metallicity tranformation based o n 95 G Cs in the Milky 
Way, M49 and M87 (Peng et al. (|2006l )). the same data 
set which is also adopted for th e colo ur-metallicity trans- 
formation by Blakeslee et al. ()2010f ). Since we do see a 
difference in the mass-metallicity relation between GCs 
in high-mass and low-mass galaxies, not only in color- 
magnitude space, it appears likely that the effect consid- 
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ered by Blakeslee et al. is in our case not the only cause 
for the environmental variation. An important step for- 
ward in this context would be to obtain accurate spectro- 
scopic metallicity measurements for an expanded sample 
(e.g., Cohen et al. (1998) and (2003)) of GCs in the 
Virgo and Fornax galaxies that define the tilt and its 
variations. 

4.2. Color Effects from Dynamical Evolution 

Having constrained the conditions under which self- 
enrichment can reproduce the observed blue tilt, it is 
worth exploring the extent to which colour changes due 
to dynamical evolution of GCs may also contribute to 
shaping the observed trends. In the context of the envi- 
ronmental variations of the tilt, this issue is worth con- 
sidering in some more detail since the dynamical evo- 
lution of star clusters is known to depend strongly on 
the external tidal field and hence varies as a function of 
environment. 

Two recent studies investigate the effects of tidal dis- 
solution of star clusters on their integrated colors (Krui- 
jssen & Lamers {2008); Anders et al. (2009)), based 
on the N-body star cluster simulations of Baumgardt 
& Makin o ((200S) . Of these two studies, Anders et 
al. (|2009( l perform a more sophisticated treatment of the 
evolution of the stellar mass function change with time, 
resulting in somewhat smaller predicted color changes. 
Both studies predict star cluster colors at a determined 
age and metallicity as a function of the ratio idyn be- 
tween current age and dissolution time. Generally, color 
changes start to set in after '^0.4 dissolution time scales. 

To determine whether these changes can be responsi- 
ble for the tilt, we estimate the amount of color change 
that has occured in the GCs of the combined VCS and 
FCS sample. We adopt a z-band mass-to-light ratio of 
M/ L^, = 1.5, independent of {g — z) color (Jordan et 
al. (|2007bD V As the fractional dynamical age, tdyn, we 
define the ratio of the current age to the dissolution time. 
In the context of the G C dissolution scenario discussed 
in Jordan et al. (|2007bf ). this can also be expressed as 
the ratio between present-day and primordial mass. Us- 
ing the present-day mass, Mpresent, and the accumulated 
mass loss, A, we find: 



ivipresent ~r ^ 

We assume A to be equal for all GCs in a given host 
galaxy, adopting A = 3 x IQ^Mq for high-mass galax- 
ies and A = 5 X IO^Mq for low-mass galaxies (see, 
e.g.. Fig. 16 of Jordan et al. ([2007b)). In the bottom 
panels of Fig. [5] we show the KMM peak positions fit- 
ted to the CMDs of GCs from the combined FCS and 
VCS high-mass and low-mass samples. We also show 
quadratic fits to the KMM peak positions. In each lumi- 
nosity bin we calculate the mean GC mass, the mean 
host galaxy magnitude and adopt A in equation |4] ac- 
cordingly, yielding an average tdyn for the GCs in each 
luminosity bin. We then calculate an estimate for the 
color change at each luminosity bin, adopting the pre- 
dictions by Kruijssen & Lamers ()2008[ ) in the V, I bands 
as a function of tdyn- We use S{V — I)/S{g — z) — 0.521 
(Peng et al. 2006) to convert to {g — z) colors. The 
results are indicated by the long-dashed magenta lines. 



It is apparent that non-negligible color changes are ex- 
pected for only the low-luminosity GCs: for > —9.0 
mag in the high-mass sample, and Mz ^ —9.5 mag for 
the low-mass sample. Moreover, the predicted strong 
changes in the red peak for the low-mass do not at 
all agree with the observational data. This comparison 
thus supports the relati vely small color changes predicted 
by Anders et al. ()2009f ): we therefore conclude that dy- 
namical evolution may be responsible for a weak color- 
magnitude trend among low-mass GCs, but certainly not 
for the strong tilt exhibited by the more massive clusters. 

We also note that the mass-loss of star clusters only 
has a very small effect on the self-enrichment estimates 
presented in the previous Section. This is illustrated by 
two example self-enrichment model calculations in the 
upper left panel of Fig.[5l The thin short dashed and thin 
solid line show the respective thick line model variations 
assuming a constant mass loss of A = 3 x lO^M© to 
have occured for all clusters. The effect on the reference 
model by Bailin & Harris is negligible, and the effect for 
the model variations that fit the data corresponds to a 
10-15% change in slope. A 0.15 dex lower pre-enrichment 
metallicity is required for the predictions that include 
mass-loss to match the predictions without mass-loss. 

5. CONCLUSIONS 

A color-magnitude relation for metal-poor GCs — the 
"blue tilt" — is now a well established observational re- 
sult. From a CMD analysis based on a sample of w 17000 
GCs detected in the ACS Virgo and Fornax Cluster Sur- 
veys, we identify a clear environmental variation in the 
slope of the GC color-magnitude relation: i.e., the slope 
of the color-magnitude relation is largest for the highest- 
mass galaxies. Nevertheless, we confirm that there are 
real, galaxy-to-galaxy differences in the amount of tilt at 
comparable galaxy luminosity and mass: e.g., M87 and 
NGC 1399 show significant blue tilts, while M49 does not 
(Mieske et al. 2006, Strader et al. 2006). 

A comparison to predictions from models of GC self- 
enrichment (Bailin & Harris (|2009[ )) suggest that this 
mechanism offers a promising explanation for the ob- 
served relation and may be able to constrain the con- 
ditions under which star formation took place in proto- 
cluster clouds. Within the context of this model, the 
observed shape of the color-magnitude relation suggests 
that protocluster clouds had density profiles somewhat 
steeper than isothermal in the highest-mass galaxies 
(/3 « 2.3). For the lower mass galaxies, a weak mass 
dependence for the structure of the protocluster clouds 
is required to match the observations (with (3 in the range 
2-5 > /? > 2.2 for masses 10^-^ < M/Mq < lO^-^^). This 
interpretation, however, is not unique. Other variants on 
this same basic picture also appear consistent with the 
observations, including models in which the clouds have 
constant density (r ~ Af^/'^) and/or relatively low star 
formation efficiencies 0.15-0.20. 

However, we point out that the significant differ- 
ences in the detailed CMDs for some galaxies of similar 
mass and morphological type pose a challenge to sin- 
gle generic mechanisms, such as self-enrichment, that at- 
tempt to explain GC color-magnitude relations. These 
differences suggest that the detailed (and stochastic) 
merger /accretion histories of individual galaxies will play 
a role shaping the appearance overall GC CMDs. 
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We conclude that further insights into the blue tilt 
phenomenon now await optical spectroscopy for large, 
representative GC samples covering the upper 3-4 mag- 
nitudes of the GC luminosity function, where the color- 
magnitude relation is most apparent. At the distances 
of the Fornax and Virgo clusters, this corresponds to ap- 
parent magnitudes of 21 < < 24.5 mag, just within 
the reach of existing 8-lOm-class telescopes. 
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Fig. 1. — Mz vs. (g — z) color-magnitude diagram (CMD) for GCs detected in the ACS Fornax Cluster Survey (FCS). From top left to 
bottom right, the plots show: (1) the GCs associated with NCC 1399; (2) all GCs in the FCS; (3) GCs belonging to high-mass galaxies; and 
(4) GCs belonging to low-mass galaxies. The vertical dashed lines in each panel mark the color expected for the blue and red populations 
based on the relation between host galaxy magnitude and peak colors presented in Peng et al. (2006). Filled circles indicate the KMM fitting 
results, corresponding to the means obtained from adopting two different initial guesses (indicated by asterisks; see §3.1 for more details). 
The number of luminosity bins is 25, except for NGC 1399, where 10 luminosity bins were adopted. The solid blue and red lines indicate 
linear least-squares fits to the KMM peak positions. The long dashed line at intermediate colors denotes the dividing point between the 
blue and red peak (i.e. the color where both blue and red peak contribute 50% to the number counts) derived from the KMM fits. The 
open triangles indicate the median color of GCs blueward and redward of the dividing line. The dotted line is a fit to the median colors. 
For both the CMD of NGC 1399 and that of the FCS high-mass sample, we also show least-squares fits to the median colors blueward 
and redward of a magnitude independent dividing line. That dividing line is given by the mean color of the long dashed line between the 
brightest magnitude bin and A/z = —8.1 mag. 
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Fig. 2. — CMDs for various GC subsamples with results from KMM overplotted. The symbols and lines are the same as in Fig.[T] (Upper 
Left Panel). The full VCS sample. (Upper Right Panel). The combined FCS and VCS sample. (Lower Left Panel). The left portion of the 
diagram shows all GCs from the high-mass FCS and VCS galaxies, while the right portion shows GCs from the low-mass FCS and VCS 
galaxies. Note the more pronounced tilt for GCs belonging to the high-mass galaxies. (Lower Right Panel): CMD showing two subsets of 
GCs in the high-mass sample, divided at a projected galactocentric distance of d = 5.5kpc into inner and outer samples. 



GC color magnitude relations 11 



NGC 1399 
1 073 GCs 



M49 
764 GCs 
N.i„=50 



M 87 

1745 GCs 



I I I I I I I I I 



T — I — I — I — I — I — I — I — p-r 
I. I 
I I ■ - 

I , i. 
I .■ I. 

I I 

r ■ -t . ' . 
. «i > .. 

■ ■ .-I . ^ "I .■ 
■I ". .. ■ ■. I . ■. 

^ r ■ . "i "-..- 




12 - 



■J- 



I 

I ■ ■ 
■J* ■ 

■ ■ 

■ I.-' 

■[.■ ■;■ 

■r * 



10 

M. 




-8 - 



-6 



■ |. . ■ ■;■ . 1 J.- . . . 



I I 

J I I I I I I I I L 



■ I 



■71 ■■■ 



■|. 



■ ■■■M:.\\r-^:/ fX.- : . 
■■.'■L.>:.'. :=. h'-^f-.^ 

■I |.- 



J I L 



1 



J I L 



1 



1 1.5 

(g-z) 



1 1.5 
(g-z) 



(g-z) 



Fig. 3. — Distribution of GCs in the CMD for three giant elliptical galaxies NGC 1399 (Fornax), M49 (Virgo), and M87 (Virgo). Results 
from KMM are overplotted for a fixed number of 15 luminosity bins. Note the absence of a blue tilt in M49, possibly due to the lack of 
intermediate-color GCs at luminosities of Mz < —10.5 mag. 
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Fig. 4. — Dependence of blue tilt on GC mass range used in the analysis. (Upper Left Panel). Variation in blue tilt slope, 7^, as a 
function of the faint magnitude cutoff, ^ut faint- In the left portion of the diagram, the filled magenta circles (green crosses) indicate the 
combined FCS and VCS sample restricted to GCs belonging to high-mass (low mass) galaxies. Open circles indicate a fit to the high-mass 
sample when disregarding the brightest GC magnitude bin at = —11.2 mag (see the bottom left panel of Fig. [2]l. In the right portion 
of the diagram, the filled green circles show the full VCS sample, while the filled magenta triangles show the full FCS sample. (Upper Right 
Panel). The mean host galaxy magnitude for GCs is plotted vs. GC magnitude, for bin size of 200 GCs. Filled green circles indicate the 
full VCS sample, while filled magenta triangles indicate the full FCS sample. The limiting magnitude between high-mass and low-mass 
galaxies is indicated by a horizontal dotted line. (Bottom Left Panel). The CMDs of GCs belonging to high-mass and low-mass galaxies. 
This is identical to the bottom left panel of Fig.[2l but with an overplotted complementary set of KMM fit results, in which we have varied 
the bright magnitude limit of the GC samples between = —12.5 mag and = —8.5 mag in steps of 0.25 mag. For each restricted 
sample, we run KMM and plot the fitted peak results. (Bottom Right Panel). This plot is based on the fits in the bottom left panel. It 
shows the slope 7z fitted to the KMM peak positions as a function of the bright magnitude cutoff, A/z cut, bright > applied to the GC samples. 
The faint limiting magnitude used in the fitting is Mz = —8.1 mag. 
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Fig. 5. — Comparison of the observed tilt (blue and red) to predicted changes of metallicity /color caused by self-enrichment (top panels) 
and dynamical evolution (bottom panels). We consider separately the GC sulssamplcs belonging to high-mass and low-mass host galaxies. 
The filled blue and red solid circles correspond to the mean GC metallicity /color in magnitude bins of 0.2 mag also shown in the b ottom left 
panel of Fig. |4] The colors (g — z) have been converted to metallicity [Z/H] using the biquadratic transformation of Peng et al. II2006I '). In 
the top left panel, the open triangles indicate the individual peak positions for NGC 1399, the squares indicate M87, and the asterisks show 
M49. (Top Panels: A Test of Self- Enrichment). The solid blue and red curves indicate quadratic fits to the KMM data points while the black 
curves are predictions from the self-enrichment model of Bail in fc Harris (2009) for various sets of model parameters, as indicated in the 
plot legend. The reference model adopted by Bailin & Harris 120091 ) is indicated by the short-dashed curve. This particular model predicts 
self-enrichment to appear at an order of magnitude larger mass than is observed. The thin short dashed and thin solid line are a variation 
of the respective thick line models, assuming a constant mass loss of A = 3 X 10^ Mq to have occured for all clusters. (Bottom Panels: A 
Test of Dynamical Evolution) . The solid curves are quadratic fits to the observed color-magnitude data points. The long dashed magenta 
curves indicate the expected color change due to dynamical evolution according to the predictions from Kruijssen & Lamers (2008). For 
the red peak at low-mass galaxies the short dashed line also shows a linear fit to the plotted data points. 



14 



Mieske et al. 



TABLE 1 

Color-magnitude trends for red and blue globular clusters in ACS Fornax and Virgo Cluster Survey galaxies, 

DETERMINED WITH KMM FITS 



Sample 


Tj.bluc 


n (Z 




A/") 




n (Z 




M") 


NGC 1399 (Mg 


= -20.9) 






-0.0S7S 


± 


0.0250 


-0.0863]* 


0.82 


± 


0.23 


-0.0256 


± 





0050 


[-0.0215]* 


0.31 


± 


0.12 


FCS high-mass (- 


22.1 < 


< -19-75 




-0.0697 


± 


0.0116 


-0.0737]** 


0.70 


± 


0.12 


-0.0363 


± 





0091 


[-0.0309]** 


0.25 


± 


0.07 


FCS low-mass (- 


19.75 < 


< -16) 




-0.01S6 


± 


0.0066 


-0.0177] 


0.28 


± 


0.10 


-1-0.0101 


± 


0.0125 


[+0.0154] 


-0.04 


± 


0.05 


FCS 








-0.0257 ± 


0.0050 


-0.0215] 


0.38 


± 


0.08 


-0.0078 


± 





0117 


[-0.0033] 


0.03 


± 


0.05 


VCS 








-0.0459 


± 


0.0048 


-0.0352] 


0.62 


± 


0.08 


-0.0223 


± 





0095 


[-0.0026] 


0.12 


± 


0.05 


FCS-I-VCS high- 


Eass (-22.1 


< Mg < ~ 


19.75) 


-0.0566 


± 


0.0076 


-0.0514] 


0.66 


± 


0.09 


-0.0198 


± 





0099 


[-0.0123] 


0.10 


± 


0.05 


FCS-I-VCS low-m 


ass (-19.75 


< Mg < - 


15.4) 


-0.0210 


± 


0.0034 


-0.0231] 


0.33 


± 


0.06 


-0.0076 


± 





0087 


[-0.0072] 


0.03 


± 


0.03 


FCS-I-VCS high- 


nass, d< 5.5 


kpc 




-0.0750 


± 


0.0120 


-0.0700] 


0.76 


± 


0.12 


-0.0301 


± 





0094 


[-0.0225] 


0.20 


± 


0.07 


FCS-I-VCS high- 


nass, d> 5.5 


kpc 




-0.0428 


± 


0.0071 


-0.0267] 


0.58 


± 


0.10 


-0.0094 


± 





0120 


[-0.0041] 


0.04 


± 


0.06 


FCS-I-VCS 








-0.0293 


± 


0.0085 


-0.0312] 


0.43 


± 


0.12 


-0.0082 


± 





0190 


[-0.0034] 


0.03 


± 


0.08 



Notes: For the samples in Column 1 (see Figs. [T] and [2] for the respective CMDs), Columns 2 and 4 give the slopes 7 between 
((? — z) and the z-band magnitudes of the blue and red GC subpopulations, determined from linear fits to the KMM peak 
positions. Errors are based on random resampling of the data points using their measured dispersion around the fit. Values 
in square brackets are the slopes determined from median colors blueward and redward of the dividing line between both 
populations. Columns 3 and 5 give the corresponding scaling relation in mass-metallicity space. The lower limit used in the 
fitting was AI^ = — 8.1 mag for all samples. 

* Slopes from median colors at magnitude independent color division: 72,biuc=— 0.0551, 72,rcd=— 0.0009 (see Fig. [T]). 
** Slopes from median colors at magnitude independent color division: 72,biuc=— 0.0331, 7z,rcd=— 0.0003 (see Fig. [T|. 



